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ABSTRACT 



We discuss identifications for 18 sources from our MAMBO 1.2mm survey 
of the region surrounding the NTT Deep Field. We have obtained accurate 
positions from Very Large Array 1.4 GHz interferometry and in a few cases IRAM 
mm interferometry, and have also made deep BVRIzJK imaging at ESO. We 
find thirteen 1.2mm sources associated with optical/near- infrared objects in the 
magnitude range K=19.0 to 22.5, while five are blank fields at K>22. We argue 
from a comparison of optical/near-infrared photometric redshifts and radio/mm 
redshift estimates that two of the thirteen optical/near- infrared objects are likely 
foreground objects distinct from the dust sources, one of them possibly lensing 
the mm source. The median redshift of the radio-identified mm sources is ~2.6 
from the radio/mm estimator, and the median optical/near-infrared photometric 
redshifts for the objects with counterparts ~2.1. This suggests that those radio- 
identified mm sources without optical/near-infrared counterparts tend to lie at 
higher redshifts than those with optical/near-infrared counterparts. Compared 
to published identifications of objects from 850/im surveys of similar depth, the 
median K and I magnitudes of our counterparts are roughly two magnitudes 
fainter and the dispersion of I-K colors is less. Real differences in the median 
redshifts, residual mis-identifications with bright objects, cosmic variance, and 
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small number statistics are likely to contribute to this significant difference, which 
also affects redshift measurement strategies. Some of the counterparts are red 
in J-K (^20%), but the contribution of such mm objects to the recently studied 
population of near-infrared selected (Js-Ks>2.3) high redshift galaxies is only of 
the order a few percent. The recovery rate of MAMBO sources by pre-selection of 
optically faint radio sources is relatively low (~25%), in contrast to some claims 
of a higher rate for SCUBA sources (~70%). In addition to this difference, the 
MAMBO sources also appear significantly fainter (~1.5 magnitudes in the I- 
band) than radio pre-selected SCUBA sources. We discuss basic properties of 
the near-infrared/(sub)mm/radio spectral energy distributions of our galaxies 
and of interferometrically identified submm sources from the literature. From a 
comparison with submm objects with CO-confirmcd spectroscopic redshifts we 
argue that roughly two thirds of the (sub)mm galaxies are at z>2.5. This fraction 
is probably larger when including sources without radio counterparts. 

Subject headings: galaxies: formation — galaxies: high-redshift — galaxies: star- 
burst — infrared: galaxies — submillimeter 

1. Introduction 

With the discovery of distant submm and mm galaxies in blank field surveys and cluster 
lens assisted surveys (see Blain et al. 2002, for a review and references), it was immediately 
recognized that this population holds important clues for the understanding of the formation 
and evolution of galaxies. A significant part of the cosmic submm background is produced 
by these objects that must be extremely luminous (Lju ~ IO^^^^^Lq) distant (z>l) infrared 
galaxies powered by intense star formation and/or powerful AGN. X-ray data argue in most 
cases against (Compton-thin) AGN and in favor of intense star formation dominating their 
luminosity (e.g., Alexander ct al. 2003). These high star formation rates (?^ilOO-1000 Mq/jt) 
and the similarity of co-moving space densities of submm sources and local ellipticals suggest 
that they are likely indicating the formation of massive spheroids. This opens a direct route to 
locating the formation of spheroids between the two extremes of an early formation similar 
to the classical 'monolithic collapse', and a late formation in hierarchical merging. More 
specifically, properties and mass functions for these high redshift objects are robust tests for 
current hierarchical models of galaxy formation (Guiderdoni et al. 1998; Kauffmann et al. 
1999; Somerville, Primack, &: Faber 2001; Baugh et al. 2003). Indeed there is evidence (e.g., 
Genzel et al. 2003; Neri et al. 2003) that these models need some modification to reproduce 
the space densities of the massive high redshift submm galaxies and their quiescent phases 
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that must exist since their star formation rates and gas content suggest a duty cycle of the 
bright (sub) mm phase well below one. 

Much of this promise will only fully come to fruition after a difficult process of identifi- 
cation and spectroscopic redshift determination. Based on optical spectroscopy, significant 
progress in the identification work was done by Chapman et al. (2003a) who presented 10 
redshifts of radio identified counterparts of SCUBA galaxies. But still, less than a dozen 
optical/near-infrared redshifts for suggested counterparts have been confirmed by CO ob- 
servations as the true redshift of the submm source (Frayer et al. 1998, 1999; Neri et al. 
2003, Greve et al. in prep, and SMM02399-0314, Kneib in prep.). Accurate positions from 
radio or mm interferometry are available for just a few dozen non-radio-preselected submm 
galaxies (e.g., Downes et al. 1999; Small et al. 2000; Eales et al. 2000; Gear et al. 2000; 
Lutz et al. 2001; Ivison et al. 2002; Ledlow et al. 2002; Webb et al. 2003a,b) and very few 
mm-selected galaxies (Bertoldi et al. 2000; Dannerbauer et al. 2002). These accurate posi- 
tions are indispensable for reliable optical/near- infrared identifications since several possible 
optical/near-infrared counterparts are usually found in the several arc second radius error 
circles of the (sub)mm surveys. This identification step is the subject of the present work. 
Of the two interferometric identification methods in use for (sub) mm sources, mm interfer- 
ometry has the advantage of directly and unambiguously locating the dust emission, but it 
is very time consuming with current instruments, with tens of hours typically needed for 
a single object in the small field of view. Because of the tight radio/far- infrared relation 
for star-forming galaxies (e.g., de Jong et al. 1985; Helou et al. 1985; Condon 1992), radio 
interferometry is the next best option, with the advantage of the large VLA primary beam 
covering one of the currently typical (sub) mm survey fields entirely. Of the brighter (sub) mm 
sources from present surveys, deep 1.4GHz VLA maps will detect all but the highest redshift 
ones (Carilli & Yun 1999, 2000a; Barger et al. 2000), and the risk of false associations of 
(sub)mm and radio sources is modest for the 1.4GHz source counts at the relevant fiux levels 
of tens of yitJy (e.g., Richards et al. 2000). 

We are building on our mm survey (Bertoldi et al. 2003, in preparation) that uses 
the MAMBO array (Kreysa et al. 1998) at the IRAM 30m telescope to cover three fields, 
the Lockman Hole, Abell 2125, and a region centered on, but larger than the NTT Deep 
Field (NDF, Arnouts et al. 1999). This paper focuses on the NDF region. Dannerbauer 
et al. (2002, hereafter Paper I) have presented results for three of the brightest NDF mm 
sources using the IRAM Plateau de Bure mm interferometer (PdBI) and BVRIK imaging. 
We now increase to 18 the number of NDF mm sources with accurate positions by using our 
VLA interferometry, and introduce additional z- and J-band optical/near-infrared imaging 
obtained since publishing Paper I. We thus identify counterparts and start assembling optical 
to radio spectral energy distributions (SEDs) of the objects constraining their nature and 
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redshift. The improved statistics is the basis for a discussion of the properties of mm galaxies 
and their relation to, or differences from, the submm population. Throughout the paper we 
adopt the cosmological parameters i/o=70kms~^ Mpc~^, 0^=0.3, Oa=0.7. 

2. Observations, Data Reduction and Results 

The radio data play a key role in the identification presented in this paper. The NTT 
deep field was observed at 1.4 GHz with the Very Large Array (VLA) in April and May, 
2001, for 15 hours in the B configuration (10 km maximum baseline). Standard wide field 
imaging mode was employed (50 MHz total bandwidth with two polarizations and 16 spectral 
channels). The source 3C 286 was used for absolute gain calibration and 1224+035 was used 
for phase and bandpass calibration. The data were also self-calibrated using sources in the 
NTT deep field itself. Images were synthesized and deconvolved using the wide field imaging 
capabilities in the AIPS task IMAGR, and the primary beam correction was applied to the 
final image. The full width at half maximum (FWHM) of the Gaussian CLEAN beam was 
7" X 5" with a PA = 0°. The rms noise on the final image is between 13 and 15 /xJy over 
the 15' xl5' region covered by MAMBO. The variation of the rms noise of the map was 
not only caused by the usual beam attenuation and UV plane coverage, but was also due to 
residual calibration errors because of a relatively large number of bright sources (2 sources of 
> 60 mJy/beam) and one spatial grouping of fainter but still relatively bright sources to the 
northeast of the center of the map. These residual calibration errors manifested themselves 
as uneven noise in a series of stripes in the final map. About 1/3 of area of the final VLA 
map was affected by this striping. We describe in § 2.1 how this was dealt with in identifying 
radio counterparts to the MAMBO 1.2mm detections. 

We use the BVRIK imaging of a field centered on the NDF that is described in Paper I. 
We have obtained additional z and J-band imaging in spring 2003. The z-band observations 
were carried out in 'service mode' on the nights of April 6 and 7, 2003 using the imaging 
spectrograph F0RS2 on UT4 of the VLT. The data were taken as a sequence of dithered 
exposures with a net integration time of 1 hour each covering 4 separate pointings. The 
pointings were chosen such that a mosaic of the images taken at the 4 pointings covers the 
whole of the K-band field centered on the NDF. The images were processed in the standard 
way but were fiat-fielded using images generated by masking out all objects with surface 
brightnesses outside 1 sigma of the background noise and then combined by taking the pixel 
median of all the frames without the images being aligned. This fiat-field frame was then 
normalized to one and then divided into each z-band image. Conditions were photomet- 
ric throughout the observations and the seeing was typically about 0.5 arc seconds. The 
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final calibration was determined through observations of the spectrophotometric standard 
LTT7987 (Hamuy et al. 1994) for the z-band filter. However, given the very sensitive red 
response of the F0RS2 array, the z-band filter was calibrated by integrating the spectral 
energy distribution of LTT7987 (Hamuy et al. 1994) and of BD+17 4708 (Oke & Gunn 
1983) which is the absolute flux calibrator for the Gunn-z filter (Schneider et al. 1983). In 
making these estimates, both the response of the CCD and the filter transmission (both of 
which are available on the ESO web pages) were considered. The average noise across the 
frame is such that the 3a detection limits in a 2" diameter aperture is zb£)+i74708=26.38. 
The absolute flux distribution of LTT7987 provides a direct calibration of the images in the 
AB- magnitude system which yields zab—'^bd+i74:708-0-^0. J-band imaging was obtained with 
SOFI at the ESO NTT between March 14 and 16 2003, covering the same 13' x 13' region 
as our Kg image of the NDF to a 3cr limiting magnitude for a 2" aperture of 24.0 in J on 
the Vega system. Data reduction followed similar steps as described in Paper 1, with the 
photometric zero point in J based on 2MASS. We used common sources to bring these new 
data to the same astrometric system as described in Paper 1, including the small correc- 
tion described in Paper I from the USNO-Al based optical/near-infrared astrometry of our 
optical/near-infrared data to the radio reference frame. All positions quoted here include 
this correction and are on the radio frame. 



2.1. Association of MAMBO sources with radio sources 

In order to obtain more accurate positions for the MAMBO sources through use of their 
radio continuum emission, we have searched the 1.4 GHz VLA data for peaks of at least 
40 //Jy (~ 3a) in the region where the MAMBO sources are located, out to radii of about 
10' from the phase center. Table 1 hsts the MAMBO sources in the catalog as of early 
2002 which have radio sources with peak positions less than 7" from the nominal MAMBO 
map position. This radius is more than half the FWHM of the MAMBO beam and will 
not exclude real associations even in case of slight systematic offsets. The table, like all 
our subsequent analysis, does not include the well-detected quasar BR1202-0725 which was 
driving the original selection of the public ESO NDF, and thus cannot be considered an 
unbiased representative of the mm galaxy population. A more complete description of the 
MAMBO data analysis will be given in Bertoldi et al. (2003, in preparation). At 5" x 7" 
resolution, our B configuration data do not put interesting constraints on the extent of the 
radio emission of the detected galaxies. Table 1 also hsts the corrected Poissonian probability 
that an association is a chance coincidence, derived using the approach of Downes et al. (1986) 
which corrects the simple Poissonian probability of an observed association for the possibility 
of associations of different nature but similar probability. In deriving these probabilities we 
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adopt the raw counts of >40//Jy peaks in the region of our VLA image covered by the 
MAMBO data, which are significantly above the true 1.4GHz source counts for a flux close 
to the detection limit. Using true radio source counts (e.g., Richards et al. 2000) would thus 
have underestimated the chance of a false association. The probabilities listed together with 
the fact that 42 sources from the original MAMBO list were searched for nearby radio peaks 
suggest that several of the associations may be chance. This statistical analysis, however, 
does not adequately consider that the excess of ~40/uJy peaks boosting the rate of chance 
coincidences is largely due to peaks in the VLA map that cluster in certain regions, due to 
uneven noise and striping because of residual calibration errors. Considering this enables 
us to identify those associations that are more hkely to be spurious by inspection of the 
distribution of >40^Jy peaks in the radio map. The quality of an associated radio source is 
indicated in Table 1, with 'uncertain' assigned to sources with an excess of similar brightness 
peaks near to the suggested radio source, and in particular with alignments of peaks due to 
striping. We also consider as uncertain radio peaks that do not remain above 40yuJy in a 
radio image where larger scale variations are partially removed by subtracting a 20" median 
smoothed version of the radio image, and two faint radio peaks near sources No. 10 and 25 
that are only above 40//Jy in this image to which we refer as the 'background-subtracted' 
image below. Of the identifications we consider 'good' below, all but source 13 are either 
above 4(7 in the original radio map, or confirmed by PdBI mm interferometry. 

Figure 1 displays the mm-radio offsets corresponding to the associations listed in Ta- 
ble 1. The dominance of real associations is clear. There were no indications for systematic 
offsets. The separations found are consistent with a typical MAMBO position error of about 

3", and confirm the radius of 7" as a sensible choice for searching for associated radio sources. 
Overall, interferometric positions are available for 18 out of 42 sources in the NDF MAMBO 
map. We consider 11 of these 18 reliable, either through 'good' radio peaks or through di- 
rect PdBI mm interferometry (Paper I). Interestingly, the PdBI data confirm the reality of 
two of the more uncertain radio sources, increasing confidence in the rest of these objects. 
The remaining MAMBO sources without radio associations could be at high redshift (with 
corresponding radio fluxes below our detection limit) , could be spurious mm sources or have 
MAMBO fluxes systematically overestimated because of the bias induced by noise in a popu- 
lation with steep number counts. Given the signal-to-noise ratios of the MAMBO detections 
all these explanations may contribute. For the purpose of this identification paper, we do 
not discuss further objects without an interferometric position. This induces a bias towards 
lower redshift and brighter objects. 
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2.2. Optical and neeir-infrcired objects neeir the radio positions 

Images of 25" x 25" size in BRzK for the MAMBO sources with interferometric posi- 
tion are shown in Figure 2. Table 2 hsts the positions and BVRIzJK photometric data of 
optical/near-infrared objects that are close to the interferometric positions. The photome- 
try is in 2" diameter apertures centered on the positions which are derived from a first step 
SExtractor (Bertin & Arnouts 1996) analysis of a noise-scaled co-addition of the BVRIzJK 
images. Wavelength-dependent morphology can lead to small offsets of the peaks in certain 
bands from these average positions (Fig. 2). We have verified the reality of objects in Table 2 
that are fainter than the 10a completeness limit and close to the detection limit by visual 
inspection of the images in the various bands. We have retained objects that are reliably 
detected individually in either B,z, or K, these images being particularly deep and including 
the extreme wavelengths first detecting blue or very red objects. For the measurements in 
the individual bands, we indicate in Table 2 cases where we consider the SExtractor aperture 
magnitudes more uncertain, since the object is only tentatively (indicated by '~') or not (in- 
dicated by '!') picked up by the eye in the visual inspection. For simplicity, we designate the 
optical/near-infrared objects with a combination of the short MAMBO number (Table 1) 
and a letter increasing with distance from the interferometric position. These letters are also 
indicated in Fig. 2. The short MAMBO numbers used in this paper are consistent with the 
NTT-MMnn numbers used by Bales et al. (2003) in their comparison of 850/xm and 1.2mm 
fiuxes of MAMBO sources. BzK color composites of the objects are shown in Fig. 3. 

We have used the publicly available hyperz photometric redshift code (Bolzonella et al. 
2000) to estimate redshifts from our seven band photometry. We have verified the hyperz 
results for our specific photometric dataset by comparison of photometric redshifts to spectro- 
scopic redshifts for 63 objects between Zspec=0.07 and 1.69 in the NTT Deep Field and its sur- 
roundings, for which spectroscopic redshifts are available from our own projects in that region 
and from ESO Science Verification data (http: //www.eso.org/science/ut2sv/NDFs_release.html). 
The dispersion (Tzspec-zphot of 0.17 is satisfactory, with few gross outliers. No trend between 
the accuracy of our photometric redshift and the color (R-K) is seen. We note however that 
hyperz does not contain an observed SED of a ultra-luminous infrared galaxy (ULIRG; which 
the SED of a mm sources may be similar), although it does allow for a range of additional 
extinction to be added to the model which may approximate a SED similar to that seen in 
ULIRGs. We notice a systematic effect at z<l, in the sense that photometric redshifts tend 
to be overestimated for sources with z<0.3 and underestimated in objects with 0.7<z<1.0. 
This is similar to what is seen in Subaru Deep Field photometric redshifts obtained with 
hyperz (Kashikawa et al. 2003, their Fig. 1). These systematic effects as well as the dispersion 
of Zspec-Zphot suggest that for individual bright objects the la errors from hyperz in Table 2 
underestimate the real errors, but at a level where the photometric redshifts are still valuable 
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constraints for the suggested counterparts and for nearby objects. 

2.3. Identification and notes on individual objects 

This section gives a case-by-case discussion of the individual interferometrically located 
MAMBO sources and identifies the suggested optical/near-infrared counterparts. We con- 
sider as possible counterparts optical/NIR objects that are less than 2 times the interfer- 
ometric positional error (Table 1) away from the interferometric position. Again, we have 
derived the corrected Poissonian probabilities for a chance coincidence computed according 
to Downes et al. (1986) and considering the search radius of twice the radio positional error. 
Table 2 lists these probabilities for the objects inside the 2a search circle, adopting the K- 
band number counts of Totani et al. (2001b). Deriving probabilities this way ignores color 
information and thus implicitly assumes a similar color distribution of counterparts and field 
population. In reality the colors differ (§ 3), i.e. we are overestimating the probability of a 
particular object being a random coinciding member of the field population. 

For convenience, salient properties of the counterparts arc collected in Table 3, full color 
information can be extracted from Table 1. Redshifts from the radio/submm spectral in- 
dex are based on the flux densities from Table 1 and the relation of Carilli & Yun (2000b). 
These redshifts are subject to the apphcability of the far-infrared SEDs used to calibrate this 
relation, and may be overestimated if cool luminous infrared galaxies (e.g. Chapman et al. 
2002b) are prominent in the (sub)mm population. The possibility of 'cool' galaxies is also 
intimately related to the role of lensing in the identification of optical/near-infrared counter- 
parts of submm sources (Blain et al. 1999; Chapman et al. 2002c) - as long as the redshift of 
the dust source is not confirmed through CO lines, certain configurations could correspond 
either to a 'cool' low z dust source or a 'warm' high z dust source with a superposed lensing 
foreground object. Therefore, for several objects, we present arguments on the potential 
role of lensing by individual foreground galaxies. These are simple plausibility arguments 
based on the observed K magnitudes and photometric redshifts, using after K-correction the 
K-band Faber- Jackson relation of Pahre et al. (1998) to estimate the velocity dispersion and 
the isothermal sphere approximation (e.g.. Peacock 1999, his eq. 4-14) to estimate the radius 
of the Einstein ring. 

MM J 120507-0748.1 (No. 03, g(ood) radio source) The faint, red counterpart 03a 

coincides with the 88;uJy radio source. A ~ 1" NW/SE double peak morphology for 3a 
may be indicated in the image but is not certain. The optical/near- IR photometric 
redshift z~2.1 is consistent with the redshift estimate based on the radio-mm spectral index 
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{zcY ~2.4). The higher estimate z=4.25 from the submm/mm ratio (Eales et al. 2003) is 
uncertain and consistent with these lower values. 

MMJ120508-0743.1 (No. 26, u(ncertain) radio source) A faint, uncertain radio 
source lies 4.6" from the nominal MAMBO position. No optical/near- infrared counterpart 
coincides with this radio source. 

MMJ120509-0740.0 (34, u) No optical/near- infrared counterpart is detected at the 
position of the faint, uncertain VLA counterpart. In the vicinity, an agglomeration of bright, 
normal galaxies is seen. The photometric redshifts of these objects (zp/jo^ RiO.4-0.5) are 
consistent with these galaxies perhaps being part of a group. At distances of >4.8" from the 
interferometric position, lensing by these individual galaxies is unlikely. Similarly, a small 
group as such is not likely to produce a strong lensing effect (e.g., Hoekstra et al. 2001). 
Eales et al. (2003) fail to significantly detect at 850/xm this object which is a Sa 1.2mm 
source in the MAMBO map and was seen at similar flux (but only 2.2a) in an independent 
short MAMBO on-off observation (Eales et al. 2003). The very low 850//m/1.2mm flux ratio 
as well as the non-detection of an optical/near-infrared counterpart may suggest a very high 
redshift. 

MMJ120510-0747.0 (10, u) Near the nominal bolometer position lies a small separation 
(1.5") galaxy pair of a fairly compact blue object (lOal) and a red (10a2) object. The 
interferometric position of the uncertain VLA source (seen only in the background-subtracted 
radio image) is about 2.4" (1.9cr) north of 10a2 and slightly more distant from lOal. It 
is thus marginally possible that one of the objects of this pair is the counterpart of the 
MAMBO source. Photometric redshifts suggest the blue object to be in the foreground 
rather than being physically associated with the red object, there arc no secondary solutions 
for the photometric redshift which would put both at the same redshift. We tentatively 
identify the source with object 10a2. The high photometric redshift, consistent with the 
radio/mm estimate, is more uncertain than its nominal error given the faintness and the need 
to decompose the two components of 10a. At the observed K-band magnitude and estimated 
photometric redshift of lOal, the estimated Einstein ring radius is only about 0.30", making 
strong lensing of 10a2 at 1.5" very unlikely. The same is true for other plausible redshifts. 
We note that lOal appears fairly compact and blue, spectroscopic observations to look for 
evidence of an AGN may be worthwhile. In case of a strong AGN its current photometric 
redshift estimate would be uncertain. 
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MMJ120516-0739.4 (36, g) This object has the strongest radio counterpart of a NDF 
1.2mm source which coincides with a relatively bright and extremely red (R-K>6.0) source. 
The discrepancy between the redshift estimate from the radio/mm spectral index (~0.8) and 
the photometric redshift (~1.7) may suggest that the radio emission is boosted by an AGN. 

MM J 1205 17-0743.1 (25, u) A faint radio counterpart is not seen in the original VLA 
image but in the background-subtracted version. It coincides with the position determined 
by mm interferometry (Dannerbauer et al. 2002, see their note added in proof). Recently, 
a faint K-band counterpart (22.5 magnitude) at the PdBI/VLA position was detected in 
ultra-deep imaging performed with ISAAC at the VLT (Lehnert et al. 2004, in preparation). 
Similarly, the new z band image indicates an extended source at the PdBI position with 
a maximum or companion to the southwest. In B, V and R very faint extended emission 
is suggested at the same position but only deeper imaging can reveal the reality of these 
structures. The photometric redshift is 2.91 in good agreement with the radio/mm redshift 
estimate and the mm/submm estimate of Eales et al. (2003). 

MMJ120519-0749.5 (01, g) Among the 'blank field' sources not showing an optical/NIR 
counterpart within 2(7 of the interferometric position, this is the blank field object with the 

strongest radio counterpart, that is also a secure radio source (72/iJy). Eales et al. (2003) 
report a A.3a 850/jm detection of this object, and infer a low 850/im/1.2mm flux ratio. This 
suggests a very high redshift which under the SED assumptions made is marginally consistent 
with the radio/mm redshift estimate. The objects a,d,c (Fig. 2) marginally detected by 
SExtractor seem to be due to noise or structure in the extended emission from the large 
foreground object and were not included in Table 2. 

MMJ120520-0738.9 (39, g) The radio counterpart is the second strongest in our sample 
and coinciding to an ERO (R-K=5.2). In the K-band, it shows weak NW/SE extensions, 
with possible variations at shorter wavelengths. The high quality photometric redshift (1-32) 
and the radio/mm estimate agree well putting the object towards the low end of the redshift 
distribution of our mm sources. Object 39f, 7.6" to the NW, has a very similar photometric 
redshift. 

MMJ120522-0745.1 (18, g) The strong radio source which is 6" SW from the nominal 
bolometer position coincides with an elongated moderate photometric redshift ERO (R- 
K=5.3), we here adopt it as counterpart of the MAMBO source noting the relatively large 
positional offset. 
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MMJ120524-0747.3 (08, u) Two faint and uncertain radio sources lie within the MAMBO 
beam. Both are blank fields in the optical/near-infrared. We adopt the /xJy radio source 
closer to the nominal MAMBO position but consider the chance of a misidentification to 
be significant, both due to the radio sources being uncertain and the difficulty of deciding 
among them (Neri et al. 2003, cf. also the case of SMM09431+4700). 

MMJ120526-0746.6 (13, g) A faint radio source from a good region of the VLA image 
coincides with an ERO (R-K=5.4). The K morphology appears distorted/curved. Both 
photometric redshift and radio/mm spectral index put the source at zai2.5. This counterpart 
has very large J-K (and z-K) color, we will discuss this below in § 3. 

MMJ120530-0741.6 (29, g) The interferometric position of the strong 85//Jy radio 
source is in a group of three objects. Object 29a coincident with the radio position is 
at I-K=3.3 a very red object in the definition of Ivison et al. (2002). 1.5" to the south-west 
there is a bluer source 29b, and 1.5" to the north-west from 29a the extremely red object 
29c (I-K>4.3). The photometric redshifts of the three objects are consistent with all being 
at z~2.4, in reasonable agreement with the radio/mm estimate of ~1.8. The three objects 
could be physically associated, with separations (for a redshift of Ri2) between about 13 and 
23kpc. We nominally identify the MAMBO source with 29a but note that all three objects 
are part of a complex structure with strongly varying colors reflecting different obscurations, 
a unique object in our sample. 

MMJ 120530-0747. 7 (07, u) A faint blue object coincides with a faint, uncertain radio 
source close to the edge of the search radius.^ It is the only candidate counterpart which is not 
detected in the K-band and the bluest in our sample. The photometric redshift of 0.35 grossly 
differs from the radio/mm estimate (~2.8), opening the possibility of the identification being 
incorrect and the object a foreground dwarf. The low mass of the relatively faint foreground 
object {Mk > —19) makes a lensing effect less hkely. The probability of ~0.1 for a chance 
coincidence (in this case based on the B magnitude and the B counts of Pozzetti et al. 1998) 
is consistent with finding such a case in our sample. 

MMJ120531-0748.1 (04, u) The identification in the optical/near-infrared is not ob- 
vious. The faint and uncertain radio source is 4.4" from the nominal MAMBO position. 



^At the chosen scaUng when making the color images, the red channel appears noisier than the blue and 
green channels. Thus, we see enhanced "red" noise in the color image. 
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behind and 1.0" from the center of an extended optical/near- infrared bright source (04a) 
with a photometric redshift estimate of 0.25. This is far below the radio/mm estimate. 
Galaxy-galaxy-lensing could hence play a significant role in MMJ120531-0748.1 — the radio 
position and the closest position on the estimated Einstein ring radius differ by less than la 
of the radio positional error. The putative lensed object must be K>;18. The morphology 
could have similarity to SMMJ04431+0210 (Small et al. 1999; Prayer et al. 2003), buth with 
a smaller separation making strong lensing more likely. 

MMJ120534-0738.3 (42, g) This is the other of the two cases where two radio sources are 
detected within our search radius. The identification of the radio counterpart seems relatively 
clear, however, since a robust brighter radio source is only 2.1" from the MAMBO position 
while the second is more distant (5.5") and is both fainter and of uncertain rehabihty. We 
adopt the first one, which is located inside a group of three faint objects (42a-c) that are 1.3 
- 2.0" (1.6 to 3a) from the radio position. An identification with 42a {zphot —l.52t.oli) or 42b 
{"^phot =2.83^024) is possible, the radio/mm redshift estimate {zcy ~2.1) is consistent with 
the photometric rcdshifts of either source. We adopt in the following an identification with 
the closest object 42a but consider an identification with b (or a blank field) possible. Source 
No. 42 appears extended in the MAMBO map, the peak fiux (Table 1) thus underestimates 
the total flux. The source is in a border region between two K-band flelds which may affect 
the quality of the K image and the photometric redshifts. 

MMJ 120539-0745. 4 (16, g) While the radio counterpart is faint, its position is con- 
firmed by the PdBI mm detection of Dannerbauer et al. (2002) which we adopt as best 
position. It is still a blank field in deeper images (K>22.7mag, Lehnert et al. 2004, in 
preparation).^ 

MMJ120545-0738.8 (40, u) A 40/xJy uncertain radio source is located 3.9" west of the 
nominal bolometer position. The faint object 40a is within 1.6a of the radio positional 
error from the radio position. Directly at the MAMBO position and 4(7 from the radio 
peak, however, a compact elliptical shaped extremely red object (40b) with a photometric 
redshift z^l.52tlit 

and a best-fitting burst-like SED is seen. An identification with this 
object may also be possible, its radio/mm index could still be consistent with the suggested 
photometric redshift range given the noise level of the VLA map. The observed radio peak 



^Thc noise peaks in the K-band seen in the color image are due to overlaping regions between two separate 
pointings in our K-band mosaic. 
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would then be unrelated. Additional MAMBO on-off observations in winter 2001/2002 could 
not conclusively settle this ambiguity but showed the strongest signal at the position of the 
ERO. Only deeper radio interferometry or mm-interferometry can clarify this ambiguity, in 
the following we adopt an identification with the ERO which is another high J-K object. 

MMJ120546-0741.5 (31, g) This source was discussed by Dannerbauer et al. (2002) 
and is a blank field at the depth of the images presented here. Lehnert et al. (2004, 
in preparation) detect a K-band counterpart (Ks=21.9mag) in VLT imaging. Eales et al. 
(2003) detect the source at 850/im and suggest a submm/mm redshift estimate consistent 
with the radio/mm estimate or with significantly higher redshifts. 

2.4. Photometric redshifts vs. radio/mm redshift estimates 

Figure 4 compares the photometric redshifts with the radio/mm redshift estimates. We 
consider this mostly a plausibility check given the uncertainties in both methods. The sources 
are often faint with relatively large photometric uncertainties in some of the band-passes (es- 
pecially the optical ones) , and the templates used in the photometric redshift estimates may 
not be fully appropriate for the likely ULIRG-hke SEDs. The radio/mm redshift estimates 
are affected by the low significance of the radio and mm detections and the internal scat- 
ter of the relation for the radio/mm redshift estimate. It is reassuring that most sources 
are consistent with the main diagonal in this diagram, with few well defined exceptions. 
Two sources (07 and 04) are at low photometric redshift but high radio/mm estimate. We 
have argued above for the possibility of foreground objects in these two cases. One object 
(36) has a relatively low radio/mm estimate which could indicate AGN contribution but 
is still consistent with scatter in this relation. The optical 'blank fields' tend to have high 
radio/mm redshift estimates. This is clearly consistent with them being faint /obscured high 
redshift sources, but mis-identifications due to the relatively uncertain radio sources may 
also contribute. 

3. Optical/neeir-infrcLred properties: very faint counterpeirts 

The optical/NIR counterparts detected for the NDF mm sources are on average faint 
even compared to published values for submm galaxies. This is illustrated in Fig. 5 which 
shows the distributions of I and K band magnitudes and limits for our sample as well as 
for three samples of submm sources where these magnitudes are available for identifications 
obtained with the help of radio or mm interferometry. For consistency with our sample, we 
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did not consider from those papers the more unrehable identifications solely based on the 
bolometer survey positions. The first of these surveys is the SCUBA 8mJy survey (Scott 
et al. 2002; Ivison et al. 2002) which, for plausible assumptions on redshift and SED of the 
objects, should be roughly similar in depth to our data. The median submm flux of the radio- 
identified 8mJy sources is Sgso^tm ~8.3mJy, while the median mm flux of our radio-identified 
sources is Si.2rram ~3.0mJy. This ratio ~2.8 of 850/xm and 1.2mm fluxes is consistent with 
z~3 ULIRG-like SEDs and is borne out in studies of well-detected robust submm and mm 
sources (e.g., Hughes et al. 1998; Ivison et al. 1998; Gear et al. 2000; Ivison et al. 2000; Lutz 
et al. 2001). We have augmented data for the second set, the SCUBA cluster lens survey 
(Small et al. 2002, and references therein), with two more objects behind one of their clusters 
detected by Cowie, Barger, & Kncib (2002) and discussed by Ledlow et al. (2002). For one 
of these wc adopt the revised identification of Ncri et al. (2003). Finally, we use magnitudes 
from the CUDSS survey, for radio identified objects in the 3h and 14h fields (Webb et al. 
2003a,b) including the object studied by Gear et al. (2000). The latter two surveys differ 
more from our survey in strategy (cluster lens) or depth (CUDSS) but may still serve as 
comparisons. 

As shown in Fig. 5, the NDF mm sources are significantly fainter both in K and in I. 
Taking the median of the observed magnitudes or limits, we find K~21.5 and I~25.5 for the 
mm objects versus K~19.6 and I~23.3 for the identifications from the 8mJy survey. The 
cluster lens and CUDSS objects have smaller statistics but show intermediate counterpart 
brightnesses (see also the recent deep SCLS imaging of Prayer et al. 2004, sometimes finding 
very faint counterparts). Ehminating all those NDF identifications based on radio sources 
that we considered more uncertain above does not change the values strongly - the median 
of the remaining objects is K~20.9 and 1^24.9, i.e. the counterparts of the MAMBO sources 
arc about 2 magnitudes fainter than those in the 8mJy 850/im survey. One major difference 
between the two populations is the absence in our sample of bright K~18 counterparts, of 
which several arc found for the 8mJy survey.^ The trivial explanation of such a difference, our 
sample just representing much fainter members of the same population, is clearly inconsistent 
with the depths and areas of the various surveys. We also do not see a trend in our sample 
for bright mm sources being brighter in K, as possibly present in the SCUBA cluster lens 
survey (Small et al. 2002). 



•^Wc note that different aperture sizes have been used in making this comparison. For example, Ivison 
et al. (2002) use a 3 arc second diameter aperture to estimate their fluxes intcad of our 2 arc second 
apertures. In sources with complex extended morphologies or nearby objects, the larger apertures could lead 
to systematically brighter magnitudes. However, it seems unlikely that this alone could lead to an offset of 
2 magnitudes. 



- 15 - 



3.1. Colors and the relation to field galetxies and the extremely red population 

Another way to represent the optical/near-infrared properties is in the form of an I-K vs. 
K magnitude-color diagram as presented for SCUBA sources in Smail et al. (2002), Ivison 
et al. (2002), and Webb et al. (2003a,b). Again, we restrict our analysis (Fig. 6) to objects 
with interferometric positions. For comparison, the figure also includes the field galaxy 
population from our NDF data and the expected properties of objects with SEDs similar to 
local ULIRGs. For these, we have picked both an object that is red in the UV/optical (Arp 
220) and a blue one (IRAS 22491-1808) from the objects studied in the UV by Trentham et al. 
(1999) and Goldader et al. (2002) and completed its SED with far-infrared data from Klaas 
et al. (2001) and large aperture optical/near-infrared data from elsewhere in the literature. 
The color-magnitude relation is then derived for these redshifted SED shapes but scaling the 
absolute flux to an observed flux of 5mJy at 1.2mm. The main results from this analysis are 
the following: First, the optical/near-infrared colors of (sub) mm galaxies scatter cover a wide 
range at a given magnitude (the 'diversity' of SCUBA galaxies, Ivison et al. 2000). Second, 
with few exceptions this range is within the envelope expected for the range of redshifted 
local ULIRG SEDs. These exceptions are from the SmJy sample which generally has a 
surprisingly large scatter of I-K colors, larger than for the mm counterparts, and includes 
some very blue objects. In our sample we flnd only one object with a limit consistent with 
I-K<3.3 and one measurement I-K=3.33, while there are about half a dozen identiflcations 
in the 8mJy sample with objects at I-K<3.3, a color well in the range of the normal fleld 
population. Third, for K>19, the average (sub)mm source is redder in I-K than an average 
field galaxy of the same magnitude. This gives some support to the suggestion that nearby 
EROs are the most hkely counterparts to SCUBA or MAMBO sources, but the dispersion in 
colors of both field galaxies and (sub) mm sources is too large to make this a robust criterion 
for identifications of individual submm sources without accurate positions. More specifically, 
the brightest (sub)mm counterparts (K<19) are statistically indistinguishable in I-K colors 
from the field population while (sub)mm counterparts are significantly redder at 19<K<21. 
It is difficult to quantify at this point how the trend continues at K>21 since observational 
limits in both K and I are significant for several of the samples involved. While the dispersion 
of I-K colors is larger for the SmJy survey counterparts than for our objects we do not see 
a significant difference in the mean color at a given magnitude. 

We compare the J-K colors which are available for our full sample with expectations for 
redshifted SEDs of ULIRGs. Fig. 7 shows the expected J-K colors as a function of redshift 
for a sample of six local ULIRGs, using their rest-frame UV to NIR SEDs assembled on 
the basis of the UV data of Trentham et al. (1999) and Goldader et al. (2002) and of large 
aperture optical/NIR photometry from the hterature. For this type of SEDs, the J-K color 
is not a unique indicator of redshift, but excursions to large J-K colors can occur for redshifts 
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above 2 and red SEDs like that of Arp220. Like Fig. 6, this diagram shows the large variety 
of colors for local ULIRGs. The locations of the MAMBO sources spread from J-K~2, 
corresponding to intrinsically bluer SEDs or lower redshifts, to J-K>3.4 which is expected 
only for red Arp220-like SEDs and high redshifts. Again, the most likely interpretation is 
that the spread in the intrinsic optical/UV SEDs of mm galaxies is large, similar to the 
spread in local ULIRGs. 

At least 5 of our 18 sources have ERO (1-K>4) counterparts. Given the long standing 
discussion about the contributions of passively evolving ellipticals and of dusty starbursts to 
this class, it is of interest to place our objects in this context. Fig. 8 shows an 1-K vs. J-K 
color-color diagram proposed by Pozzetti & Mannucci (2000) to address this issue. Based on 
the dissimilarity between a 4000A break and a reddened more smooth starburst continuum, 
this diagram is expected to separate the two categories for redshifts in the range of roughly 
1 to 2. With the possible exception of MMJ12507-0748.1 (No. 03) the proposed mm and 
submm counterparts avoid the region of z~l-2 evolved populations, which might be found 
in case of false identifications of mm sources with such objects. It is not possible to draw the 
inverse conclusion from this diagram, i.e. to conclude that we are observing starbursts at 
z~l-2. The region populated by the (sub) mm sources is indicative of smooth spectra which 
are not very specific on redshift, or of spectra with a break between J and K. This caveat 
is particularly relevant since some (sub)mm objects are significantly fainter than the z~l-2 
EROs for which this method was developed. 

Recently, Franx et al. (2003) presented a technique to select optical-break galaxies at 
a redshift z>2 with the near-infrared filters and K^. van Dokkum et al. (2003) reported 
spectroscopic redshifts of a handful of such (J5-K5)>2.3 pre-selected galaxies, five of six 
indeed lying at redshifts higher than z>2 (z=2.4-3.5). Individual similar objects have been 
found in other surveys (e.g., Dickinson et al. 2000; Maihara et al. 2001; Totani et al. 2001a; 
Im et al. 2002, sometimes called 'HEROs'). We have probed for an overlap between this 
population and the MAMBO identifications. Franx et al. (2003) report a surface density of 
3 ± 0.8 arcmin"^ for objects with Js-Ks>2.3 and <22.5. In our sample from a survey 
area of roughly 100 square arcmin we suggest 3 identifications with J-K>2.3, K<22 galaxies. 
Given several more J>24 limits, a few more might be found in deeper images to K=22.5 
among the objects with current K>22 hmits, and among the MAMBO sources without 
radio identification that are not discussed in this paper. It is thus fair to assume that the 
contribution of mm galaxies in the few mJy range to the population of luminous red galaxies 
as observed by Franx et al. (2003) is of the order one to a few percent. This small overlap is 
consistent with submm/mm followup for the four HEROs of Totani et al. (2001a), reporting 
very weak sub/mm signals (Coppin et al. 2002, and Andreani et al. in prep.). Also, only 
one of the six Js-Ks>2.3 objects with spectroscopic redshifts reported by van Dokkum et al. 
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(2003) has a spectrum consistent with that of a dusty starburst. On the other hand, the 
fraction of red J-K sources (J5-K5>2.3 and <22.5) in the mm population is significantly 
higher. Of our MAMBO sources with radio counterparts, at least 3 out of 18 (~20%) are 
identified with red J-K sources. It may be higher, but the lack of deeper J-band data prevents 
us from estimating the colors of the faintest of the K-band counterparts. There might be 
an evolutionary connection beyond this small overlap, i.e. a contribution of mm galaxies in 
quiescent or post-active phases to the population of luminous red galaxies/ HEROs. 

3.2. Dependence on Radio Properties: Radio-Pre-Selection and Detection 

Limits 

One of the major issues in our developing understanding of the nature of (sub) mm emit- 
ting galaxies at high redshift is how does the method of selection influence the properties of 
the sources? For example, targeted submm follow-up of optically faint radio sources (OFRS: 
SiAGHz > 40yuJy, I>25; Chapman et al. 2001) has been used as an efficient method to detect 
and study a significant 'radio pre-selected' fraction of the (sub)mm population. This is a 
particularly important comparison for our study. Because we have limited ourselves to a 
investigation of only those mm sources in the NDF with interferometric identifications, we 
might expect a radio-preselected sample to have similar properties. Indeed, the strength of 
the conclusions from such studies obviously depends on the amount of overlap with the pop- 
ulation detected in unbiased 'blind' surveys (i.e., those not dependent of radio pre-selection), 
and on the nature of the biases induced by the OFRS selection. Both Barger et al. (2000) 
and Chapman et al. (2001), in comparing the submm detection efficiency of radio sources 
with that of blank field surveys, report a recovery rate of about 70% of bright (>5-6mJy) 
submm sources when targeting Si^ghz ^ 40//Jy radio sources with very faint near-infrared 
or optical counterparts (HK' >20.5 or I>25, respectively). In contrast, only 10% of the 
refined 8mJy sample are reported to have an Siaghz > 40/iJy, I>25 counterpart (Ivison et 
al. 2002). Using the same criteria, we recover 10 of 42 (about 25%) of the original NDF 
MAMBO fist. 

While the recovery fraction of optically faint radio identified sources presented here is 
similar to the low SmJy recovery rate of Ivison ct al. (2002) at first glance, a clear difference 
exists in the nature of the sources not matching the Siaghz > 40;uJy, I>25 criteria. We 
find 7 to 8 (about 20%) optically bright (I<25) radio sources, while Ivison et al. (2002) find 
12 of 30 optically bright sources (40%) in their refined sample with Siaghz > 40/iJy radio 
counterparts. This is another manifestation of the clear difference in typical counterpart 
brightnesses between the NDF MAMBO sources and the SmJy survey reported above. Ir- 
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respective of these differences, a large recovery rate by the radio pre-selection technique for 
faint optical sources is not supported by both 'blind' surveys. 

Part of the difference between the suggested 70-75% recovery rates and the lower num- 
bers found in the blind surveys may be simply due to reference integral counts of the complete 
submm population and due to small areas that we used in these surveys. The counts adopted 
by Barger ct al. (2000) and Chapman et al. (2001) are more towards the low end of the spread 
of the count determinations summarized in Blain et al. (2002) and the number of sources in 
these field meeting their criteria was relatively small (for example, in Chapman et al. (2001), 
the recovery rate was determined based on 11 radio sources with submm detections and 2 
sources without radio counterparts). In any case, a large recovery rate of the optically faint 
radio pre-selection is not ensured, and the selection effects are complex (see also Chapman 
et al. 2002a). While the radio detection will overall favor low redshift (CariUi & Yun 1999), 
optical faintness will prefer among those objects the higher redshifts ones. This is indeed 
suggested in our radio-identified sample where the mean redshift (approximated by the ra- 
dio/mm indicator) is < z >— 2.8 ± 0.2 for the I>25 objects compared to < z >— 2.4 ± 0.2 
overall. 

While a comparison of the recovery rate of radio pre-selected sources with faint optical 
counterparts is interesting and instructive, perhaps a more robust comparison can be made 
with Chapman et al. (2003c). In Chapman et al. (2003c), the radio pre-selected sample 
is defined without imposing constraints on magnitudes of the optical/near-infrared coun- 
terparts to the radio sources, now finding noticeable numbers of optically brighter (1<25) 
radio-identified submm sources in contrast to the earlier radio-preselected studies. Over- 
all, they recover about 70% of all submm sources which is formally higher than both the 
MAMBO source recovery rate for our survey and that of Ivison et al. (2002, which have 
recovery rates of ~40% and ~60% respectively). Similar to the differences we have found 
between our results and those of the 8 mJy SCUBA survey of Ivison et al. (2002), we also 
find a substantial difference between our results and the sources meeting our selection crite- 
ria of Ssso^tm ~8mJy and Siaghz > 40y[xJy in Chapman et al. (2003c). For the 10 sources in 
Chapman et al. (2003c) meeting these criteria, we find a median I-band magnitude of ~24 
with only 20% of the sources having I-band magnitudes fainter that 1=25 (2 out of 10). Such 
a result is in stark contrast to our results with a median I-band magnitude of ~25.5 and 
~80% of sources having I>25. 

It is important to note that the differences in recovery rates and counterpart properties 
are not a direct consequence of the applied radio fiux limits combined with strong trends 
with flux. Making a higher flux cut (40// Jy as in our sample) in the radio detections of Ivison 
et al. (2002) only removes a small number of sources from their survey (3 or 4 depending 
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on how to treat a binary radio source, leaving a recovery rate of still ~50%). Similarly, the 

majority of sources in the radio preselected samples of Barger et al. (2000), Chapman et al. 
(2001), and Chapman et al. (2003c) have Si^ghz > SOyuJy (well above their detection limit), 
changing the lower radio flux to our value thus has little effect . 

4. Optical to radio SED constraints 

We discuss the near-infrared (rest frame optical) to radio spectral energy distributions 
of our sources using an updated version of the diagram presented in Paper I, combining 
a spectral index between the K band and the sub mm and the well known radio/submm 
spectral index (Fig. 9). To be able to compare MAMBO and SCUBA sources, we estimate 
850^m fluxes for those MAMBO sources where they have not been measured, by scaling 
the 1.2mm fluxes by a factor of 2.5 which is consistent with the median value of published 
ratios. This implicitly assumes that there is no major difference between the submm and 
mm populations. In Fig. 9, the MAMBO objects are also on average below the SCUBA 
ones, consistent with their optical/near- infrared faintness derived in the previous section. 

Given the opposite trends in K corrections for the (sub)mm flux on one side and the 
radio and optical/NIR fluxes on the other side, redshifting a dusty galaxy SED will move 
the object from top left to lower right in this diagram. Such a trend is indeed indicated 
in the upper panel. The lower panel verifies and helps to visualize this expectation by 
plotting the loci for the redshifted SEDs of six local ULIRGs for which UV to radio SEDs 
are available. These clearly provide an indication of the scatter in the properties of the 
local population. In particular, significant offsets to the left are possible for objects with 
strong ACN, such as IRAS 19254-7245 with its reddened Type 2 AGN (Mirabel et al. 1991) 
and strong radio emission. It will be interesting to probe for AGN activity for the objects 
in the left (Q;(850;Um,1.4GHz)<0.65) part of the diagram, indeed some of the objects there 
show evidence for AGN activity (LE850.12, Ivison et al. (2002), J02399-0134, Soucafl et 
al. (1999)). For the radio/submm spectral index on the horizontal axis the dispersion for 
a larger sample is of the order 0.16 (CariUi & Yun 2000a), for the NIR/submm index of 
similar magnitude given the scatter arising from adopting different ULIRG SEDs (Fig. 3 of 
Paper I.). A similar scatter is observed in the high redshift (sub)mm population, shown 
by labeling in the top panel objects with CO-confirmed redshifts with the redshift value. 
Four of these six objects may also serve as an indication of a locus of z~2.5 objects that 
is independent of the uncertain applicability of local templates. The average locus of these 
z=2.4 to 2.8 (mean z=2.6) objects is a(850^m,1.4GHz)=0.76 and a(2.2^m,850//m)=-0.99. 
Better statistics from on-going CO follow-up observations providing confirmed redshifts is 
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highly desirable to improve this calibration of Fig. 9. Although, the CO confirmed redshift 
may introduce an upward bias for this point if the successful optical redshift measurements 
on which they depend on correspond to brighter counterparts on average. Less bias should 
be present in the horizontal (radio/sub mm) direction. We find about two thirds of the 
MAMBO/SCUBA objects towards the 'high redshift' (z>2.5) direction to the right of this 
point (43/69 at Q;(850/xm,1.4GHz)>0.76). This is an underestimate, since the plot shows 
almost exclusively radio detected objects, thus missing roughly one-third to a half of the 
parent samples including the potentially highest redshift objects. Altogether, this indicates 
that the median redshift of 2.4 suggested for the bright SCUBA population from optical 
redshifts (Chapman et al. 2003a) is an underestimate of the median redshift of the complete 
population. 

A conclusion drawn in Paper I from this type of diagram is that MAMBO and SCUBA 
sources are at high redshift and/or are at least as obscured as local ULIRGs. A technical 
factor affecting this reasoning is due to the unknown spatial structure of the objects. If the 
objects were very large and considering the effects of cosmological surface brightness dim- 
ming, the observed small structures (Fig. 2) may just refiect the central high surface bright- 
ness regions of more luminous extended objects that cannot be fully retrieved at realistic 
noise levels by just increasing aperture size. This would cause a tendency to underestimate 
the total near-infrared flux. There is evidence for some objects of that type in the (sub)mm 
population from near-IR imaging (e.g., Lockman 850.1, Lutz et al. 2001), and radio map- 
ping (Ivison et al. 2002) finds some objects resolved by a 1.4" beam. Examples of extended 
and complex morphologies are also seen in HST images of radio-preselected submm galaxies 
(Chapman et al. 2003b). For many other well studied (sub)mm sources the detected regions 
fit into 2" apertures but this may be deceiving effects of surface brightness dimming if they 
are at very high redshift. Some estimate of the magnitude of the effect can be gained from 
comparison to the local ULIRG population. Half light radii have been estimated at several 
wavelengths, with median results of the order 5.3 kpc in U (Surace & Sanders 2000), 3 kpc in 
I (Zheng et al. 1999) and about 1.5 kpc in H and K (Cohna et al. 2001; Tacconi et al. 2002). 
At redshifts around 3 an aperture of 2" diameter spans about 15 kpc. Half light radii of ~4 
kpc in the rest frame V band (observed K band) would then mean no significant loss of light 
outside our aperture. This is only true as long as the structural analogy to local ULIRGs 
holds, of course, and a major fraction of light could be lost for objects more than a factor 
of 2 larger, especially if at z>4. Such large objects would be placed too low in Fig. 9 by our 
measurements. This implies that the displacement of many SCUBA and MAMBO objects 
down-wards from the location of local ULIRGs could include an effect of large spatial extent 
as well as obscuration of the rest frame optical part of the SED. Testing such a scenario 
would require extremely deep K imaging. Size constraints at other wavelengths (optical. 
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mm, radio) will be useful as well, but for this question they will have to be interpreted with 
some caution considering the strong trends for local infrared galaxies, where half light radii 
decrease as wavelength increases. 

5. What causes the difference between SCUBA and MAMBO populations? 

The difference in brightness of about two magnitudes between typical MAMBO and 
SCUBA counterparts is a surprising result, given the selection of these objects at fairly 
similar wavelengths and on a part of the rest frame spectral energy distribution that does not 
show prominent breaks or features. We discuss in the following several potential contributors 
to this difference. 

Given the faintness of the targets in comparison to instrumental sensitivities at all 
wavelengths, the identification process may still contain errors. For the NDF MAMBO 
sources, some of the near 3a radio sources are clearly uncertain and may lead to uncertain 
identifications. We have argued above, however, that eliminating the uncertain sources as 
classified in § 2.1 does not change the magnitude distribution of the MAMBO counterparts 
significantly. For the SCUBA sources similar problems may be present at the limit of the 
radio data. In addition, there are by now examples where intensified identification efforts 
have rejected previously adopted bright counterparts and lead to a significant upward revision 
of the proposed counterpart magnitudes, the most notable being HDF850.1 (Hughes et al. 
1998; Richards 1999; Downes et al. 1999; Dunlop et al. 2002), SMMJ00266+1708 (Frayer et 
al. 2000) and SMMJ09431+4700 (Neri et al. 2003). Taking away a few of the brightest or 
bluest 8mJy counterparts from Fig. 5 does not fully remove the difference to the MAMBO 
counterparts, though. Noting all these uncertainties we do not see any compelling evidence 
for systematic mis-identifications. Identification uncertainties are also related to the presence 
of 'complex' objects with several components, e.g., one of the mm galaxies discussed by 
Bertoldi et al. (2000). Ivison et al. (2002) emphasize the role of such objects in the 8mJy 
survey. Again, the example of SMMJ09431+4700 (Neri et al. 2003) where the submm and 
CO source H7 is 4" from the object H6 with a very similar (optical) redshift indicates that 
such grouping or clustering exists. In such cases, the (sub) mm source could relate to just 
one or to all of the components. We find only three such objects in our sample. For No. 10 
we have preferred one component because of photometric redshift arguments. For objects 
29 and 42 there is more uncertainty on whether to identify the source with one or several 
components, but even the sum of the components will not be a 'bright' counterpart. 

Optical/near-infrared magnitude differences could also reflect a significant difference 
in the redshift distribution of the two populations. The effect in the optical/near-infrared 
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will obviously depend on the intrinsic SEDs, for ULIRG-like SED shapes and keeping the 
(sub)mm flux fixed a difference in median redshifts of the order Az=l-2 may be necessary 
(Fig. 6). Chapman et al. (2003a) suggest a median redshift of 2.4 for bright submm galaxies 
from optical spectroscopy. Some of these optical redshifts still require confirmation by CO 
detections, and it is natural to assume that the optical redshifts are still biased towards the 
optically bright and low redshift end of the population. This assumption is supported by 
SED arguments (§ 4). Even at redshifts around or somewhat above 3, however, both SCUBA 
and MAMBO observations effectively sample the long wavelength side of the rest frame far- 
infared SED peak, unless the intrinsic dust temperatures were extremely cold. Evidence for 
both wavelengths sampling the long wavelength side is indeed present in many of the best 
studied SCUBA sources (e.g., Ivison et al. 1998; Downes et al. 1999; Gear et al. 2000; Ivison 
et al. 2000; Lutz et al. 2001). At such redshifts, it will be difficult to construct a redshift 
distribution that fully accounts for the observed differences. Redshift estimates for the 
SCUBA population from the radio/submm spectral index (Carilli & Yun 2000a), and from 
estimates trying to more fully include the rest-frame far-infrared part of the SED under the 
assumption of SEDs similar to local ULIRGs (Yun & CariUi 2002; Wiklind 2003; Aretxaga et 
al. 2003) indicate median redshifts in the range 2.5 to 3.5 with the possibility of a significant 
high redshift tail (Wiklind 2003; Aretxaga et al. 2003). However, it is important to note 
that such analysis suffer from possible template mismatches between high and low redshift 
sources and using an unrealistically narrow range of SED types and properties when making 
redshift estimates. For example, Aretxaga et al. (2003) may perhaps put undue weight on 
single templates by matching one template to each source to constrain the possible redshift 
range. This means that low redshift SEDs they have in their analysis must be good analogues 
for high redshift sources. In addition, a substantial number of higher redshift z>4 sources in 
the MAMBO population is also suggested by the 850//m/1.2mm ratios presented by Eales et 
al. (2003), albeit with large uncertainties due to the technical difficulties of the method. For 
the objects overlapping with our study (Nos. 1, 3, 16, 25, 31) their 850/Ltm/1.2mm redshift 
estimates can be compared to the detection in various optical pass-bands - if they are at high 
redshift, these objects should be undetected dropouts in optical bands. With the possible 
exception of No. 3 (best submm/mm estimate z~4.25 but a marginal detection in B) the 
results are consistent, and even for this object considering the uncertainties. In summary, 
while there may be a high redshift tail of the (sub) mm population preferentially sampled 
by MAMBO, we find it difficult to explain the magnitude differences by a bulk redshift 
difference. This view is consistent with Fig. 9 which does not only show the 'diagonal' offset 
between the two sets expected in such a case. 

Because of the degeneracy of redshift and dust temperature in interpretations of the 
submm/radio spectral energy distributions (e.g. Blain et al. 2003), similar arguments can be 
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made with respect to a possible separation of the MAMBO and SCUBA populations by dust 
temperature, for example due to an evolutionary difference between MAMBO and SCUBA 
objects. One should note, however, the need for significant obscuration in the rest frame 
optical when invoking 'cool' rest- frame FIR SEDs (e.g.. Paper I). 

The results shown here indicate that there are real differences in the counterpart bright- 
nesses and perhaps redshifts of the SCUBA and MAMBO populations. On the sole basis 
of the SED arguments presented it is difficult to robustly weight the contributions of real 
differences in redshifts and spectral energy distributions of the mm and submm selected pop- 
ulations to this difference, in comparison to caveats due to small sample statistics, cosmic 
variance, and identification difliculties. Significant progress in settling these issues will come 
from tests of proposed identifications with available optical redshifts through CO measure- 
ments (e.g., Neri et al. 2003), and from future direct submm/mm redshift measurements for 
the fainter part of the population, most likely through wide band searches for CO emission. 

6. Conclusions 

We have discussed optical/near-infrared identifications for those 18 of 42 sources in our 
MAMBO 1.2mm map of the NTT Deep Field region for which interferometric positions are 
available through a VLA 1.4GHz map and in three cases IRAM PdBI mm interferometry. 
In addition to being the basis for identifications, our deep BVRIzJK imaging allows the 
derivation of optical/near- infrared photometric redshifts for the counterparts and for nearby 
objects. In comparison with radio/submm redshift estimates, these photometric redshifts 
suggest that two of the optical/near-infrared sources close to interferometric positions are in 
the foreground, in one case likely lensing the background mm source. One strongly lenscd 
object in this sample is consistent with expectations for the (sub)mm population (Blain et 
al. 1999; Chapman et al. 2002c). This leaves us with eleven detections of counterparts at 
magnitudes of K=19 to 22.5, and seven limits or blank fields with most hmits at K>22. The 
I-K and J-K colors of the counterparts are consistent with redshifted SEDs similar to local 
ultra-luminous infrared galaxies, and likely with a similarly large spread of the rest frame 
UV/optical SED properties. The counterparts of mm sources contribute to the recently 
discussed population of J-K>2.3, K<22.5 high redshift galaxies, but only at the few percent 
level for the current mm survey depths. At least ~20% of the MAMBO sources with radio 
counterparts have J-K>2.3 and K<22.5. 

The counterparts to NDF mm sources are on average about 2 magnitudes fainter than 
counterparts presented for the 8mJy 850/im survey which is of similar depth and are similarly 
faint compared to radio pre-selected submm sources. Remaining mis-identifications, redshift 
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or temperature differences between tfie two populations, and small number statistics/cosmic 
variance all may contribute to this difference, at levels that are hard to quantify from cur- 
rently available data. Our result reinforces the view that direct (e.g., wide band CO) spectro- 
scopic redshifts may be necessary for a substantial fraction of the (sub)mm population which 
is very faint in the optical/near-infrared. From a comparison of near-infrared/submm/radio 
spectral indices with those of submm sources with CO-confirmed redshifts we suggest that 
the fraction of (sub)mm galaxies at z>2.5 is about two thirds for the interferometrically 
located ones and larger when adding the radio-undetected part of the population. 

We would like to thank the staffs of ESO, IRAM, and VLA for their support with the 
data acquisition and reduction. We thank Andrew Baker for helpful discussions. We also 
acknowledge the contributions of the referee in clarifying a number of important points. This 
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by the Jet Propulsion Laboratory, California Institute of Technology, under contract with 
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Table 1. NDF MAMBO sources and associated 1.4GHz radio sources 
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Comment 



(12) 



PdBI pos. 



PdBI pos. 



PdBI pos. 



Note. — Col. (1) — MAMBO source. 
Col. (2) — Short number on the original MAMBO NDF source list. 

Col. (3)-(4) — ,72000 coordinates of the source in the MAMBO data. 

Col. (5) — MAMBO 1.2mm peak flux density. *: The source MMJ12U534-0738.3 appears to be extended and may be multiple or smeared. For consistency with 
the other sources, we list its peak flux which is effectively a lower limit to the true flux. The source was included because of a reliable detection in total flux and 
a 1.4GHz counterpart. 

Col. (6)-{7) — J2000 coordinates of the associated radio source. 

Col. (8) — VLA 1.4GHz peak flux density. The average rms of the radio map is 13 /iJy. See also col. (11) for possible effects of striping residuals caused by 
bright sources contained in the radio map (see § 2). 

Col. (9) — Separation between MAMBO position and VLA position. 
Col. (1(3) — Probability that the association is a chance coincidence. 

Col. (11) — Quality flag for good (g) and uncertain (u) radio sources. To declare the quality of a radio source as good two criteria has to fulfilled: (1) source 
does not disappear in the 'background-subtracted' radio map, and (2) no eye-catching clustering of >40/iJy peaks in the adjacent region. 
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Table 2. Optical/near-infrared objects close to the interferometric positions 



Source RA DEC Sep B V R I z J K ^phot P 

(J2000) (J2000) arcsec mag mag mag mag mag mag mag 
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Table 2 — Continued 



Source RA DEC Sep B V R, I z J K ^-phot P 

(J2000) (J2000) arcscc mag mag mag mag mag mag mag 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 
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Table 2 — Continued 



Source RA DEC Sep B V R, I z J K ^-phot P 

(J2000) (J2000) arcscc mag mag mag mag mag mag mag 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 
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Table 2 — Continued 



Source RA DEC Sep B V R, I z J K ^-phot P 

(J2000) (J2000) arcscc mag mag mag mag mag mag mag 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 



Table 2 — Continued 



Source RA DEC Sep B V R, I z J K ^-phot P 

(J2000) (J2000) arcscc mag mag mag mag mag mag mag 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 
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Table 2 — Continued 



Source RA DEC Sep B V R, I z J K ^-phot P 

(J2000) (J2000) arcscc mag mag mag mag mag mag mag 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 
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Table 2 — Continued 



Source RA DEC Sep B V R, I z J K ^-phot P 

(J2000) (J2000) arcscc mag mag mag mag mag mag mag 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 



Table 2 — Continued 



Source RA DEC Sep B V R, I z J K ^-phot P 

(J2000) (J2000) arcscc mag mag mag mag mag mag mag 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 
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Table 2 — Continued 



Source RA DEC Sep B V R, I z J K ^-phot P 

(J2000) (J2000) arcscc mag mag mag mag mag mag mag 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 
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Table 2 — Continued 



Source RA DEC Sep B V R, I z J K ^-phot P 

(J2000) (J2000) arcscc mag mag mag mag mag mag mag 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 
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Table 2 — Continued 



Source RA DEC Sep B V R, I z J K ^-phot P 

(J2000) (J2000) arcscc mag mag mag mag mag mag mag 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 
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Table 2 — Continued 



Source R.A DEC Sep B V R I z J K ^phot P 

{J2000) (J2000) arcsec mag mag mag mag mag mag mag 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 
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4 


,2 


>27.4 


>26.5 


>26.2 


!25.4 


>26.4 


>24.0 


21.04 


4 


.21+°- 


16a 


12 


;()5 


;39, 


,61 


-07; 


;45; 


;27, 


,9 


2 


,3 


25.54 


25.30 


24.85 


23.63 


24.32 


»23.8 


21.11 





,93tg: 


16b 


12 


()5 


;39 


,35 


-07; 


;45; 


;28, 


,7 


2 


,4 


25.47 


25.27 


24.78 


23.88 


24.95 


22.21 


20.50 


1 




16c 


12 


()5 


39 


,54 


-07; 


;45; 


;24, 


,3 


2 


,9 


25.20 


24.58 


23,59 


23.06 


23.59 


21.97 


20.45 





,46t°; 


16d 


12 


;()5 


;39, 


,19 


-07; 


:45 


:22 


,7 


6 


,0 


26.21 


K26.3 


a26,0 


24.62 


25.16 


= 23.5 


=21.9 


1 


,i3tS; 


16e 


12; 


:05; 


:39, 


,93 


-07; 


:45; 


:25, 


,0 


7, 


,1 


26.05 


!26.4 


fB25.4 


S325.0 


25.65 


>24.0 


>22.0 





,58lS; 


40a 


12; 


;06; 


;46, 


,69 


-07; 


;38; 


;60, 


,6 


2, 


,1 


26.33 


ra26.1 


«25.9 


!26.0 


!26.1 


!23.7 


!21.7 


1 


,6itS; 


40b 


12; 


;06; 


;46, 


,98 


-07; 


;38; 


;61, 


,8 


3, 


,9 


>27.4 


>26.5 


>26.2 


!26.3 


K326.8 


22.19 


19.57 


1 




40c 


12; 


;06; 


;45, 


,66 


-07; 


;38; 


;48, 


,4 


4, 


,8 


f»26.7 


>26.5 


«25.2 


!26.2 


26.40 


>24.0 


>22.0 







40d 


12; 


;06; 


;46, 


,69 


-07; 


;38; 


;69, 


,1 


6, 


,6 


25.24 


24.63 


23.92 


23.60 


23.99 


!23.2 


21.41 


3 




31a 


12; 


;06; 


;46, 


,69 


-07; 


;41 


;30, 


,3 


4, 


,3 


25.80 


25.04 


«24.8 


Ri23.8 


26.12 


!23.6 


ra21.8 


2 


,87t»„; 


31b 


12; 


;06; 


;46, 


,94 


-07; 


;41 


;33, 


,3 


5, 


,3 


25.52 


25.70 


25.33 


!«24.7 


26.08 


K323.1 


>22.0 


1 


,35lg; 


31c 


12; 


;05; 


;46, 


,66 


-07; 


;41; 


;39, 


,6 


5, 


,3 


25.72 


25.23 


24.48 


23.30 


23.76 


22.29 


20.37 


1 




31d 


12; 


;06; 


;46, 


,82 


-07; 


;41 


;38, 


,9 


5, 


,7 


«26.8 


>26.5 


25.01 


23.82 


24.68 


K323.0 


>22.0 


0, 


,70±0; 



0.060 
0.082 



Note. — Col. (1) — Optical/ncar-infrarcd source identification. Suggested counterparts to mm sources are marked in bold. 
Col. (2)-(3) — J2000 coordinates of optical/near-infrarcd source. 

Col. (4) — Separation between optical position and best interferometric position. 

Col. (5)-(ll) — Source magnitudes in a 2" aperture. All magnitudes are on the Vega system except for the z magnitudes which are on the BD 
system (see text). 

Col. (12) — Optical/near-infrared photometric redshift. 

Col. (13) — Probability of a chance coincidence with an object of the given K magnitude (07a: B) at the given separation. 

f : magnitude contaminated by nearby foreground galaxy. 
*: magnitude contaminated by nearby foreground star. 

*: 03d is not detected in B,z and K but retained due to credible R and I detections. 

o: magnitude contaminated by lOal. 

• : Lehnert et al. (2004, in preparation). 
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Table 3. Main properties of counterparts to MAMBO sources 



(1) 



(2) 



Q 

(3) 



Sl.2™„ 
mjy 
(4) 



Si, 4 

(s) 



mag 
(6) 



(7) 



"^■phot 
(8) 



(9) 



MMJ120607-0748. 


1 


03 


g 


4, 


.6 


88 


MMJ120608-0743. 


1 


26 


u 


2, 


.7 


42 


MMJ120609-0740. 


,0 


34 


u 


3, 


.3 


47 


MMJ1206 10-0747. 


,0 


10 


u 


3, 


.1 


41 


MMJ120616-0739. 


,4 


36 


g 


2, 


.2 


469 


MMJ120517-0743. 


,1 


25 


g 


4, 


.3 


40 


MMJ120619-0749. 


,5 


01 


g 


5, 


.2 


72 


MMJ120520-0738. 


,9 


39 


g 


2, 


,7 


336 


MMJ120522-0745. 


1 


18 


g 


2 


,0 


90 


MMJ12()r)24-()747, 


3 


08 


u 


2 


,0 


47 


MMJ12()526-()74(i 


6 


13 


g 


2 


,6 


41 


MMJ12()53()-0741 


6 


2!) 


g 


2 


,3 


85 


MMJ12()53()-0747. 


7 


07 


u 


3 


.0 


40 


MMJ12()531-0748. 


1 


04 


u 


2 


,7 


53 


MMJ12()534-0738. 


,3 


42 


g 


2 


2 


61 


MMJ120539-0745 


4 


16 


g 


3 


,4 




MMJ120545-0738. 


,8 


40 


u 


3, 


.7 


40 


MMJ120646-0741. 


,5 


31 


g 


6, 


.6 


42 



21.48 
>22.0 
>22.0 
20.58 
20.11 
22.50 
>22.0 
19.04 
19.61 
>22.0 
20.44 
20.68 
>22.0 
>18.0 
20.94 
>22.7 
19.57 
21.90 



2 64"*' -^'^^ 
^■°*-0.95 
„ 70 + 1.49 

2-»2_ j gg 

„ „„ + 0.42 
U-»J_o.34 

-1.08 
2.74+' 



1.05: 

1.671 
2.2lj 
2.62+ 



-0.87 
c + O.SO 
-0.39 
7+0.79 
-0.68 
1+1.25 
-0.79 
,+ 1.64 
-0.92 
, „„ + O.SS 
l'**^-0.63 
, 00+1.77 

„ ov+1.32 
^"^'-0.86 

2 05+1" 
^■"■'-0.79 

^•■"'-0.84 

3 05+1-'"' 
■^■"f -1.04 

3.86l?J« 



9 n+0.26 



2 49+0 °l 

^- '^-0.20 
2.91+" 



1+0.12 
-0.14 



1.32 
1.25 



+0.04 
0.01 
+0.06 
0.06 



K+0.21 
-0.59 

2.26iS:;« 

(0.35tS:f,) 

1 t.,+0.24 
^■°^-0.34 



2.26"f 



,+1.24 
''-0.25 



4 25+~ 

(see note) 



3.66 



+4.26 
2.53 
> 8.96 



_+6.14 
"-1.93 



Note. — Col. (1) — MAMBO source. 
Col. (2) — Short number on the original MAMBO NDF source list. 

Col. (3) — Quality of the interferometric identification, g(ood) stands for objects with a good quality VLA 

source or a PdBI identification. 

Col. (4) — MAMBO 1.2mm peak fiux density. 

Col. (5) — VLA 1.4GHz peak fiux density. 

Col. (6) — K magnitude (Vega), see Table 2 for other bands. For MMJ120619-0749.5 and MMJ120631-0748.1, 
wc give K-band limits as we assume that the foreground objects detected at the radio position are not linked to 
the dust emission. 

Col. (7) — Redshift estimate from radio/mm spectral index (Carilli & Yun 2000b). Errors reflect Ict uncertainties 
in alpha due to measurement errors and due to the intrinsic scatter of the relation, added in quadrature. 
Col. (8) — Photometric redshift from optical/near-infrared photometry. Values in brackets refer to objects - 
MMJ120619-0749.5 and MMJ120531-0748.1 - that are likely in the foreground. 

Col. (9) — Redshift estimate from the 1200/im/850/um fiux ratio as derived by Eales et al. (2003). No. 34 was 
observed but the flux ratio could not be flt by the templates used at any redshift. 
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Fig. 1. — Histogram of separations between mm positions from the MAMBO map and VLA 
positions of associated radio sources. Light and dark grey indicate the good and uncertain 
sources from Table 1. The dashed histogram indicates the expected number of unrelated 
background radio peaks. 

Fig. 2. — BRzK images (from left to right) for the fields of the MAMBO sources with radio 
associations. The 25" x 25" images are centered on the nominal MAMBO positions hsted 

in Tabic 1 and oriented such that north is at the top and cast to the left. A small cross 
indicates the best intcrfcromctric position (PdBI if available, VLA otherwise) and its la 
error. The right panel indicates all > 40yuJy radio peaks by crosses, and the 7" radius used 
to search for radio peaks possibly associated with the MAMBO sources. Small letters label 
optical/near-infrared sources near to the best interferometric position (see also Table 2). 

Fig. 3. — BzK color composites of the NDF MAMBO sources with interferometric positions. 
Large 7" radius circles are drawn around the nominal MAMBO positions. Unlike in Fig. 2, 
we now show both the VLA 1.4GHz peaks (indicated by squares) and position of the PdBI 
mm interferometric sources (indicated by small circles). 

Fig. 4. — Comparison of optical/near-infrared photometric redshifts with redshift estimates 
from the radio/mm spectral index. Blank fields without photometric redshifts are shown 
with an arbitrary lower limit of 0.5 for the optical/near- infrared photometric redshift, and 
without the error bars for the radio/mm estimate. For two objects shown by squares we 
believe that the measured optical/NIR photometric redshift refers to a foreground object, 
the true redshift of the submm source being higher (see text). 

Fig. 5. — Distribution of K and I magnitudes for the interferometrically located NDF 1.2mm 
objects and for three samples of interferometrically located SCUBA 850/xm galaxies. The 
limit symbol at K>22, I> 25.5 in the NDF panel stands for 5 objects. In the panel for the 
8m Jy survey, we indicate the regions occupied in this type of diagram by the class O/I/II 
objects in the definition used by Ivison et al. (2002). In the top two panels thick symbols 
reflect identifications described as good/robust. Magnitudes for the cluster lens sources have 
been corrected for the magnifications stated in the original papers. 

Fig. 6. — I-K vs. K magnitude-color diagram for the MAMBO identifications (large filled 
circles). Small dots show the field galaxy population from our NDF data. Interferometri- 
cally located SCUBA galaxies are represented by squares for the 8m Jy sample objects from 
Ivison et al. (2002), stars for cluster lens survey objects (corrected for magnification) from 
Small et al. (2002), and triangles for the CUDSS objects of Webb et al. (2003a,b). Filled 
diamonds connected by continuous lines show the expected colors/magnitudes for redshift 
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1,2,... for objects with the SED shapes of Arp220 and IRAS 22491-1808 but scaled to ob- 
served Si,2mm=5mJy. Average properties in magnitude bins are shown with error bars for the 
(sub)mm population (open diamonds) including all objects shown, and for the field popula- 
tion (crosses). Blank field sources are not included, the figure is an incomplete representation 
in its right part. 

Fig. 7. — J-K color as a function of redshift for the template SED shapes of a number 
of local ultra-luminous infrared galaxies. Horizontal over-plotted lines show measurements 
(continuous) or lower limits (dashed) for the J-K color our MAMBO sources. 

Fig. 8. — I-K vs. J-K diagram, with the thick dashed lines showing the regions populated by 
z~l-2 passive eUipticals and by dusty starbursts according to Pozzetti & Mannucci (2000). 
Triangles indicate NDF MAMBO sources. Diamonds indicate SCUBA sources from the 
literature with UK photometry and one A2125 MAMBO source from Bertoldi et al. (2000). 
The thin vertical dashed line indicates the J-K>2.3 criterion of Franx et al. (2003). 

Fig. 9. — Top panel: Diagnostic diagram for the near-infrared to radio SEDs of (sub)mm 
sources, combining the radio/submm spectral index with a spectral index between near- 
infrared and submm. Thick squares/limits represent the NDF MAMBO source discussed 

in this paper, thin asterisks/limits mostly SCUBA sources from the samples described in 
the text, and a few individual well studied (sub)mm objects. Objects with CO-confirmcd 
redshifts are labeled with their redshift. A typical observational error is indicated in the 
lower left. The lower panel shows for comparison the loci corresponding to the redshiftcd 
SEDs of a number of local ULIRGs with available UV to radio SEDs. Note that IRAS 
19254-7245 hosts a powerful AGN. 
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